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We review the perturbative QCD formalism for exclusive heavy meson 
decays, concentrating on the three-scale factorization theorem for nonleptonic 
processes. The formalism is then extended to the radiative decay B — > K*^, 
which occurs through penguin diagrams. It is observed that the contributions 
from the operators other than the penguin one b — > S7 are not negligible. 
From the best fit to the experimental data of the branching ratio B(B — > 
K*^), we extract the B meson wave function, which possesses a sharp peak 
in the region with a small momentum fraction. 



1 talk presented at the Fourth Workshop on Particle Physics Phenomenology, Kaohsi- 
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1. Introduction 



Recently, perturbative QCD (PQCD) has been proposed to be an alter- 
native approach to the study of heavy hadron decays [1], ||, Q , which com- 
plements the heavy quark effective theory |4]] and the Bauer- Stech-Wirbel 
(BSW) method |J. In this talk we review the PQCD analysis of exclu- 
sive B meson decays, concentrating on the three-scale factorization theorem 
for nonleptonic processes ||. The PQCD formalism is then applied to the 
radiative decay B — > K*^, which involves the contributions from penguin 
diagrams. According to the factorization theorem, a heavy meson decay rate 
is expressed as the convolution of a hard subamplitude with meson wave 
functions. The former is calculable in usual perturbation theory, while the 
latter must be extracted from experimental data or derived by using non- 
perturbative methods, such as QCD sum rules. Since the K* meson wave 
function has been known from the sum rule analyses |7] , the B — > K *7 decay 
is an ideal process, from which the unknown B meson wave function can be 
determined. With the B meson wave function obtained here, we are able to 
make predictions for other B meson decay modes. 

The operator for exclusive semileptonic heavy meson decays, such as B — > 
D^lv, is 

^V^Q. - rilcrfb , (1) 

where Gf is the Fermi coupling constant and V c b the Cabibbo-Kabayashi- 
Maskawa (CKM) matrix element. The formalism for semileptonic decays is 
simpler, and part of that for nonleptonic decays. Hence, we shall explain 
only the factorization of nonleptonic heavy meson decays below, which occur 
through the Hamiltonian, 

H=^V t3 V: i {q l q k ){q ] q t ), (2) 

(qq) = 57^(1 — 75)5 being the V — A current. Hard gluon corrections cause 
an operator mixing, and their renormalizat ion-group summation leads to the 
effective Hamiltonian, 

tfeff = ^WSIdCAOOiM + c 2 (/i)0 2 (//)] , (3) 
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with the four-fermion operators 0\ = (qiqk)(QjQi) an d 2 = (qjqk)(<ilQi)- The 
Wilson coefficients Ci >2 , organizing the large logarithms from the hard gluon 
corrections to all orders, describe the evolution from the W boson mass My/ 
to a lower scale \i with the matching conditions Ci(M w ) = 1 and c 2 (M w ) = 0. 

The most widely adopted approach to exclusive nonleptonic heavy meson 
decays is the BSW model ||, in which the factorization hypothesis on the 
matrix elements of the operators 0\ 2 is assumed. In this model decay rates 
are expressed in terms of various hadronic transition form factors. Employ- 
ing the Fierz transformation, the coefficient of the form factors corresponding 
to external W boson emissions is a± = c\ + C2/N, and that corresponding 
to internal W boson emissions is a 2 = c 2 + C\/N, N being the number of 
colors. The form factors may be related to each other by heavy quark sym- 
metry, and parametrized by different ansatz. Nonfactorizable contributions, 
which can not be expressed in terms of hadronic transition form factors, and 
nonspectator contributions from W boson exchanges are neglected. 

Though the BSW model is simple and gives predictions in fair agreement 
with experimental data, it encounters several difficulties. It has been known 
that the large N limit of ai.2, i.e., the choice a\ = C\(M C ) « 1.26 and 
a 2 = c 2 (M c ) « —0.52, with M c the c quark mass, explains the data of charm 
decays ||. However, the same large N limit of ai = ci(M&) ~ 1.12 and 
a 2 = c 2 (Mb) ~ —0.26, Mb being the b quark mass, does not apply to the 
bottom case. That is, the different mechanism between charm and bottom 
decays can not be understood in the BSW approach. Moreover, it has been 
difficult to explain the two ratios associated with the B — > J/vpK^ decays 



simultaneously using the BSW model, where B(B — > J/ipK) is the branching 
ratio of the decay B —>■ J / ipK. This controversy was found to be attributed 
to the neglect of the nonfactorizable internal iy-emission amplitudes, which 
are in fact of the same order as the factorizable ones [|, [Hfl . 

The above difficulties have been overcome by a modified PQCD formal- 
ism, the three-scale factorization theorem 0. In the next section we shall 
summarize the basic ideas of this theorem. 
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B(B -> J/ijK*) 
B(B -»• J/ipK) ' 



Rl 



BjB -> J/jjKj) 
B(B -> J/ijK*) ' 



(4) 
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2. Three-scale factorization theorem 



Nonleptonic heavy meson decays involve three scales: the W boson mass 
Mw, at which the matching conditions of the effective Hamiltonian to the 
original Hamiltonian are defined, the typical scale t of a hard subamplitude, 
which reflects the dynamics of heavy meson decays, and the factorization 
scale 1/6, with 6 the conjugate variable of parton transverse momenta. The 
dynamics below 1/6 is regarded as being completely nonperturbative, and 
parametrized into a meson wave fucntion <p(x), x being the momentum frac- 
tion. PQCD for the scale above 1/6 is reliable, and radiative corrections 
produce two types of large logarithms ln(M^/t) and \n(tb). The former are 
summed to give the evolution from My/ down to t described by the Wilson 
coefficients c(t). While the latter are summed to give the evolution from t to 
1/6. 

There exist also double logarithms ln 2 (P6) from the overlap of collinear 
and soft divergences, P being the dominant light-cone component of a meson 
momentum. The resummation of these double logarithms leads to a Sudakov 
form factor exp[— s(P, b)], which suppresses the long-distance contributions 
in the large 6 region, and improves the applicability of PQCD around the 
energy scale of few GeV. The 6 quark mass scale is located in the range of 
applicability. This is the motivation we develop the PQCD formalism for 
heavy hadron decays. For the detailed derivation of the Sudakov factor for 
heavy meson decays, refer to |l], 0. With all the large logarithms organized, 
the remaining finite contributions are absorbed into a hard 6 quark decay 
subamplitude H(t). Because of Sudakov suppression, the perturbative ex- 
pansion of H in the coupling constant a s makes sense. 

Therefore, a three-scale factorization formula possesses the typical ex- 
pression, 



c(t) <g> H(t) <g> <p(x) <g> exp 



-s(P,b)-2 f ^ lq (a s (fl)) 

Jl/b jJL 



(5) 



where the exponential containing the quark anomalous dimension ^ q de- 
scribes the evolution from t to 1/6 mentioned above. As the hard scale t runs 
to below Mb and M c , the constructive and destructive interferences involved 
in bottom and charm decays, respectively, appear naturally. Furthermore, 
not only factorizable, but nonfactorizable and nonspectator contributions 
can be evaluated in a systematic way. With the inclusion of nonfactorizable 
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contributions, we find that a± j2 restore their original role of the Wilson coef- 
ficients, instead of being treated as the BSW free parameters. The branching 
ratios of the various decay modes B — > D^Tc(p) and the ratios R and Rl 
associated with the B — > J/tjjK^ decays can all be well explained by the 
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We review the evaluation of the factorizable, nonfactorizable and non- 
spectator contributions to exclusive nonleptonic decays B — > D^7r(p) in 
the three-scale PQCD factorization theorem. The factorizable external W- 
emission amplitudes define the B — > D^*> transition form factors i = +, 
— , V, A±, A 2 , and A3, through the hadronic matrix elements, 



(D{P 2 )\V^B{P 1 )) = y]M B M D [£+{ri){v x + v 2 )» + - 



(D*(P 2 )\V^\B(P 1 )) = isJM B M D .£ v (T ] )e^ e * V2aVip ? 



-U(rj)<? • ihvfl . (6) 

Pi (P2), ^ffi {M D (*)) and t>i (^2) are the momentum, the mass, and the 
velocity of the B (D^) meson, satisfying the relation Pi = MgVi (P; 
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M D (*)V 2 ). e* is the polarization vector of the £)* meson. The velocity transfer 
V\ ■ v 2 in two-body nonleptonic decays takes the maximal value 77 = (1 + 
r 2 ) / (2r) with r = M D (*)/Mb- In the infinite mass limit of and M D {*) the 
form factors £j obey the relations 

t + = iv = u = u 3 = e, e- = a 2 = o, (7) 



where £ is the so-called Isgur-Wise (IW) function |L2 . 

Equation (^) is the standard definition of the form factors In our 
approach, however, the factorization formulas for will contain the Wilson 
coefficients as the extra convoluiton factors. £j include the contributions from 
the hadronic matrix element of 0\ and from the color-suppressed matrix 
element of 2 . Therefore, their factorization formulas involve the Wilson 
coefficient a\ = c\ + c 2 /N. The form factors for the internal TV-emission 
diagrams include the factorizable contributions from the matrix elements of 
2 and from the color-suppressed matrix element of Oi, which then contain 
the Wilson coefficient a 2 = c 2 + Ci/N. The form factors for the IV-exchange 
diagrams include the factorizable contributions from the matrix elements of 
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O2 and from the color-suppressed matrix element of 0\. Hence, they also 
contain the Wilson coefficient 

The nonfactorizable contributions to the B — > D^*'tt decays are color- 
suppressed, which can not be expressed in terms of hadronic form factors. 
A nonfactorizable diagram associated with external W emissions involves a 
hard exchanged gluon, which, for example, attaches the spectator quark of 
the B meson and a valence quark of the pion. The resultant amplitude de- 
pends on the Wilson coefficient C2/N. A nonfactorizable diagram associated 
with internal W emissions contains a hard gluon, which attaches the specta- 
tor quark of the B meson and a valence quark of the D meson. The resultant 
amplitude depends on the Wilson coefficient c\/N. Similarly, a nonfactoriz- 
able amplitude associated with W exchanges contains the Wilson coefficient 
Ci/N. For the explicit expressions of the above form factors and amplitudes, 
refer to UTT . 



3. The B -> K*^ decay 

We now apply the three-scale factorization theorem to the radiative decay 
B — > fT*7 |L3| |, whose effective Hamiltonian is 
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H cS = -^*E^)°iM • ( 8 ) 
V 2 i=i 



It will be found that only the operators 



2 = (sc)(cb) , (9) 



e 



7 = —M b sa^{\ + 75)^ , (10) 



in 



o 8 = - JLm 6 <^(i + i 5 )t%g; u , (11) 

with o^ v = (i/2)[7 M , 7"], are important. The \x dependence of the corre- 



sponding Wilson coefficients 02(11), c-j(^l) and cs(n) can be quoted from [[14 

We assign Pi = (M B /V2)(1, 1, Or) and P 2 = (M B /y/2)(l, r 2 , T ) with 
r = Mk*/Mb, Mk* being the K* meson mass, as the B and K* meson 
momenta, respectively, in the rest frame of the B meson. Assume that the 
light valence quark in the B (K*) meson carries the momentum k\ (hi)- k\ 
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has a minus component fcf, giving the momentum fraction x B = k^ / P{ ', and 
small transverse components k#T- k 2 has a large plus component k 2 , giving 
xr — k 2 I ' P 2 ', and small \s-kt- The photon momentum is then P 3 = Pi — P 2 , 
whose nonvanishing component is only P 3 ~ . 
The decay amplitude is expressed as 

M = e 7 • e* K .Mi + ie w+ .t^,M 2 , (12) 

with e 7 and e#* the polarization vectors of the photon and of the K* meson, 
respectively, and 

M m = M 1(2)2 + M 1(2)7 + M 1{2)8 , (13) 

where the terms on the right-hand side come from the operators 2 , O7, 
and Os, respectively. From Eq. fll2"|) , it is obvious that only the K* mesons 
with transverse polarizations are produced in the decay. The factorization 
formulas for the amplitudes My are, in terms of the overall factor, 

r<°> = %-v: s v th c F Ml , (14) 

listed below. The B — > K*^ decay rate is then given by 

_ 1 -r 2 



(|Mx| 2 + \M 2 \ 2 ) . (15) 



32ttM b 

The amplitude from 2 is written as 

. 4 Z" 1 pl—x pi fl/A 

M 12 = P - 1 - dx dy dxsdxK \ bdb(p B (x B )(f)K*{xK) 
3 jo jo jo jo 

xa s (t 2 )c 2 (t 2 ) exp[-S(x B , x K , t 2 , b, b)} 
x [(1 — r 2 + 2txk + 2xb)v ~ ( rx K + 3:tb)(1 — x)] 
(1 — r)(l — t 2 )xkx 



xy(l - r 2 )x K M 2 B - M 2 ' 



H 2 (Ab,J\B 2 \b), (16) 



with 



A z = x K x B M 2 B 



,2 _ _ ^2 y ,,2 , M c 



5 2 = x K x B M B - -^(l-r z )x K M B + 



1 — x x(l — x) 



max(A,J\B 2 \,l/b) , (17) 
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where the integrations over x and y are associated with the charm loop, from 
which the photon and the hard exchanged gluon are radiated. The variable 
ipx), conjugate to the parton transverse momentum k B j> (^at), represents 
the transverse extent of the B (K*) meson. t 2 is the characteristic scale of 
the hard subamplitude, 



as mentioned before. 

There are two diagrams associated with the operator 7 , where the hard 
gluon attaches the two quarks in the B meson and in the K* meson sepa- 
rately. The corresponding amplitude is 



H 2 (Ab, J\Bj\b) 



K (Ab) - K (J\Bj\b) , B 2 2 >0, 



K (Ab) - i^H£\J\Bj\b) , Bl < , (18) 



Myj 




Jo Jo 



x rHi a \Ab K ,^\b B ,J\Bj\b K )F 7 (t 7a ) 

+ [l + r+(l-2r)x K ]Hi b \Ab B ,C 7 b B ,C 7 b K )F 7 (t 7b )} , (19) 



with 



B 2 = (x B -r 2 )M 2 B , C 2 = x K M 2 




t 7b = max(A, C 7 , l/b B , l/b K ) . 



(20) 



The function F 7 denotes the product 



F 7 (t) =a s (t)c 7 (t) exp[-S(x B ,x K ,t,b B ,b K )} . 



(21) 



The hard functions 



Hi a) (Ab K , J\Bj\b B , J\Bj\b K ) = 



K (Ab k )h(J\Bj\b B , y/\Bf\b K ) , B 2 > , 



K (Ab k )ti(J\Bj\b B , J\Bj\b K ) , B 2 < , 



(22) 
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with 

h = 0(b B -b K )K o (y/\Bj\b B )I o (J\Bj\b K ) + (b B ^b K ) , 



ti 

and 



7T 

l 2 



H?\Ab B , C 7 b B , C 7 b K ) = K (Ab B )h(C 7 b B , C 7 b K ) , (24) 



9{b B - b K )H^\J\Bj\b B )J Q (J\Bj\b K ) + (6 B <- b K 



(23) 



come from the two diagrams. 

Four diagrams are associated with the operator 8 , where the photon 
is radiated by each quark in the B and K* mesons. The corresponding 
amplitude is 

M 18 = -r (0) - / dx B dx K / b B db B b K db K (j)B(x B )(j)K*(xK) 
o Jo Jo 

x |(1 - r 2 + x B )(r^ + x B )^ a) (A6^, £ 8 6 B , B 8 b K )F 8 (t 8 a) 

+ [(2 - 3r)a;^ - x B + r(l - x K )(rx K - 2rx B + 3x B )] 

xH^{Ab B ,C 8 b B ,C 8 b K )F 8 {t 8b ) 
+ (1 + r)(l - r 2 )[(l + r)2; B - rx K ] 

xHi c \J\AV\b K , D 8 b B , D 8 b K )F 8 (t 8c ) 
_[(! _ r 2)((i _ r 2)( 2 _ XK ) + (1 + 3r)(2a* - x B )) 
+2r 2 x K (x K - 

xH { 8 d \J\A^\b B , E 8 b B , E 8 b K )F 8 (t 8d )\ (25) 



with 



A' 2 = (1 - r 2 )(^ - ^)M 2 , = (1 - r 2 + x B )M% , 

C 2 8 = (1 - x x )Ml , 7J 2 = (1 - r 2 )x B M 2 B , El = (1 - r 2 )x x Ml , 

t 8a = max(A, B 8 , l/b B , l/b K ) , 

t 8b = max(A, C 8 , l/b B , l/b K ) , 

t 8c = max( v ^, D 8 , l/b B , l/b K ) , 



t 8d = m a x(J\A'%E 8: l/b B: l/b K ) . (26) 



9 



The funciton F 8 denotes the product 

F 8 (t) = a s (t)c 8 (t) exp[-S(x B ,x K ,t,b B ,b K )] . 

The hard functions 

Hi a) (Ab K , B 8 b B , B 8 b K ) = H?\Ab K , B 8 b B , B 8 b K ) , 
Hi b) (Ab B , C 8 b B , C 8 b K ) = K (Ab B )h'(C 8 b B , C 8 b K ) 
Hjf\J\A^\b K ,D 8 b B ,D 8 b K ) = 

A' 2 > , 



K (^\A'i\b K )h(D 8 b B ,D 8 b K ) 
i\HS\J\A»\b K )h{D 8 b B ,D 8 b K ) , A' 2 <0 , 

H<£\J\A^\b B ,E 8 b B ,E 8 b K ) = 



K (J\A»\b B )ti(E 8 b B ,E 8 b K ) , A' 2 > , 



qH^(^/\A^\b B )h'(E 8 b B , E 8 b K ) , A' 2 < 

are derived from the four diagrams. 

The amplitudes with the structure ie^-e^e*^* are 



M- 
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r o / / d y dx B dx K / bdb4> B (x B )4>K*( x K) 
6 Jo Jo Jo Jo 

xa s (t 2 )c 2 (t 2 ) exp[-S(x B , x K , t 2 , b, b)} 
x "((1 - r)(l - r 2 ) + 2r 2 ^ + 2x B ) y 
-(r(l + r)x K + (3-r)x B )(l-x)} 



x 



(1 — t 2 )xkx 



M 27 
M 28 



xy(l -r 2 )x K M 2 B - M c 2 
-M 17 , 

l/A 



1 /-l /-l/A 

r (0) - / dx B dx K / b B db B b K db K (p B (x B )(f) K *(x K ) 
6 Jo Jo 

x |(1 - r 2 + a; B )(ra* + x B )^ a) (^, £ 8 & B , B 8 b K )F 8 (t 8a ) 

+ [(2 - 3r)x^ - - r(l - x K )(rx K - 2ra;B + 3x B )] 
xHi\Ab B ,C 8 b B ,C 8 b K )F 8 (t 8b ) 
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+ (1 — r)(l - r 2 )[(l + r)x B - rx K ] 

xHi c \J\A»\b K , D 8 b B , D 8 b K )F(t 8c ) 
_[(! _ r 2)((i _ r 2) (2 + ^) _ (i _ 3r ) XB ) 

-2r 2 x K (^ - xb)) 
rid) 



xH { 8 a> (J\A' 2 \b B ,E 8 b B ,E 8 b K )F(t 8 



(31) 



It is obvious that the above factorization formulas bear the features of Eq. (|5|). 

The exponentials exp(— S) appearing in My are the complete Sudakov 
factor, with the exponent 



S = s(x B P^~,b B ) + 2 — j q (a 8 (n)) 

Jl/bn u 



/b B H 



+s(x K P+, b K ) + s((l - x K )P+, b K ) + 2 f -H 7 (a a {p)). (32) 

Jl/b K H 



l/b K fl 

The wave functions <fr B and <pK* satisfy the normalization, 



(33) 



with the decay constant fi, i = B and K* . The wave funciton for the K* 
meson with transverse polarizations has been derived using QCD sum rules 
H, given by 

fa. = V6f K *x(l - a;) [0.7 - (1 - 2x) 2 ] . (34) 



As to the B meson wave functions, we investigate two models |15| |16 

N B x(l - x) 2 



4> ( b\x) 

4>f\x) 



M B + C B (l-x) 



N' B Jx(l — x) exp 



1 / xM B 

2 V u 



(35) 
(36) 



with the normalization constants N B and N' B , and the shape parameters C B 
and ijj. 



4. Results 
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In the evaluation of the various form factors and amplitudes, we adopt 
G F = 1.16639 x 1(T 5 GeV~ 2 , the decay constants f B = 200 MeV and f K * = 
220 MeV, the CKM matrix elements \V t *V tb \ = 0.04, the masses M c = 1.5 
GeV, M B = 5.28 GeV and M K * = 0.892 GeV, and the B° meson lifetime 
r B o = 1.53 ps [11]. It turns that no matter what value of the shape parameter 
C B is chosen, the model (p B in Eq. fl3"5| ) with a flat profile leads to results 
much smaller than the CLEO data of the branching ratio BiB — > K*^) = 
(4.2 ± 0.8 ± 0.6) x 10~ 5 [17|: the maximal preiction from 4> B \ corrsponding 
to the shape parameter C B = —M B , is about 1.1 x 10~ 5 . The model <fi B 
in Eq. ( |3^| ) with a sharp peak at small x then serves a possible candidate. 
It is indeed found that as u = 0.425 GeV, a prediction 4.19 x 10~ 5 for 
the branching ratio is reached, which is equal to the central value of the 
experimental data. If varying the shape parameter to u = 0.42 GeV and to 
u = 0.43 GeV, we obtain the branching ratios 4.36 x 10" 5 and 4.08 x 10~ 5 , 
respectively. It indicates that the allowed range of u is wide due to the large 
uncertainties of the data. 

The detailed contribution from each amplitude Mij is listed in the follow- 
ing table in the unit of 10 -6 GeV~ 2 : 



M 12 /r 


M 17 /r 


M 18 /r 


-6.20 - 20.38i 


11.10 - 198.08i 


-1.94-0.29i 


M 22 /r 


M 27 /r 


M 28 /r 


-4.84- 16.17* 


-11.10 + 198.08i 


2.13 + 0.36i 



It is observed that the contributions from the operator 2 are comparable 
with those from the penguin operator O-j. The amplitudes associated with 
0$ are less important. The other operators give contributions smaller than 
those from 8 by an order of magnitude, which are thus not included here. 
Furthermore, the imaginary parts of Mum from the timelike interal quarks 
play an essential role for the explanation of the data. 

In this work we have determined the B meson wave function, 



b B (x) = 1.81359 v /x(l - x) exp 



xM 



B 



0.425 GeV 



(37) 



from the best fit to the experimental data of B(B — > K*y). We stress that, 



however, Eq. (pf\ ) is not conclusive because of the large allowed range of the 
shape parameter uj. A more precise B meson wave function can be obtained 
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by considering a global fit to the data of various decay modes, including 
B — > D^n(p) ||18|| . Once the B meson wave fucntion is fixed, we shall 
employ it in the evaluation of the nonleptonic charmless decays. 

This work was supported by the National Science Council of Republic of 
China under the Grant No. NSC-88-2112-M006-013. 
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